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^ ■ We present X-ray and optical spectroscopic observations of twelve galaxy groups 

Qh and clusters identified within the XMM Large-Scale Structure (LSS) survey. Groups 

q , and clusters are selected as extended X-ray sources from a 3.5 deg 2 XMM image 

5^ ■ mosaic above a flux limit 8 x 10~ 15 ergs -1 cm~ 2 in the [0.5-2] keV energy band. 

| Deep BVRI images and multi-object spectroscopy confirm each source as a galaxy 

C3 concentration located within the redshift interval 0.29 < z < 0.56. We combine line-of- 

' sight velocity dispersions with the X-ray properties of each structure computed from 

. £h ! a two-dimensional surface brightness model and a single temperature fit to the XMM 

. spectral data. The resulting distribution of X-ray luminosity, temperature and velocity 

!h ' dispersion indicate that the XMM-LSS survey is detecting low-mass clusters and 

galaxy groups to redshifts z < 0.6. Confirmed systems display little or no evidence for 
X-ray luminosity evolution at a given X-ray temperature compared to lower redshift 
X-ray group and cluster samples. A more complete understanding of these trends will 
be possible with the compilation of a statistically complete sample of galaxy groups 
and clusters anticipated within the continuing XMM-LSS survey. 

Key words: X-rays: galaxies: clusters; Cosmology: large-scale structure of the uni- 
verse; Surveys 



ABSTRACT 



1 INTRODUCTION 

Surveys of distant galaxy clusters map the distribution in 
the universe of large amplitude density fluctuations, and 
so constrain key cosmological parameters and permit sec- 
ondary studies to determine how X-ray gas and galaxy evo- 

* Based upon observations performed at Paranal (70.A-0283), 
Las Campanas and CTIO observatories and on observations ob- 
tained with XMM-Newton, an ESA science mission with instru- 
ments and contributions directly funded by ESA Member States 
and NASA 
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lution proceeds as a function of environment. Wide area 
X-ray surveys are well placed to compile statistically well- 
defined samples of distant galaxy clusters because a) both 
source confusion and the X-ray background are low com- 
pared to optical searches, b) computation of the selection 
function and volume sampled is straightforward and c) the 
selection of extended X-ray emitting sources is sensitive to 
the signature of hot gas contained within massive, gravita- 
tionally bound structures. 

A number of systematic X-ray studies have extended 
both the maximum redshift (i.e. the most luminous galaxy 
clusters) and the minimum luminosity (i.e the least mas- 
sive structures) to which X-ray clusters can be identified. A 
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comprehensive review is provided by Rosati, Borgani & Nor- 
man (2002). The principal aim of such surveys for distant, 
X-ray emitting clusters is to determine their space density 
evolution as a function of redshift and to constrain the com- 
bination of the root mean square mass density fluctuations 
on 8/1 ~ Mpc scales, as-, and the overall matter density of 
the Universe expressed as fraction of the closure density, Om 
(e.g. Borgani et al. 2001; Schuecker et al. 2001, Allen et al. 
2003). In addition to the study of global cosmological pa- 
rameters, galaxy clusters provide examples of dense cosmic 
environments in which it is possible to study the evolution of 
the hot, X-ray emitting gas (e.g. Ettori et al. 2004; Lumb et 
al. 2004) and to determine the nature of the cluster galaxy 
populations and the physical processes underlying observed 
trends in galaxy evolution (e.g. Yee et al. 1996; Dressier et 
al. 1999; Andreon et al. 2004a). 

Extending our current knowledge of low luminosity (i.e. 
low mass) X-ray clusters represents an important challenge 
for the present generation of X-ray surveys performed with 
the XMM and Chandra facilities. The local (z < 0.25) X-ray 
Luminosity Function (XLF) for galaxy clusters is currently 
determined to X-ray luminosities Lx ~ 10 42 ergs -1 in the 
[0.5,2] keV energy band (Henry et al. 1992; Rosati et al. 
1998; Ledlow et al. 1999; Boehringer et al. 2002). However, 
our understanding of such systems at redshifts 0.25 < z < 
0.8 is largely restricted to luminosities Lx > 10 43 ' 5 ergs -1 
(Henry et al. 1992; Burke et al. 1997; Rosati et al. 1998; 
Vikhlinin et al. 1998; Ebeling et al. 2001). 

Low luminosity (10 42 < L x /ergs -1 < 10 43 ) X-ray 
clusters correspond to low mass clusters and larger galaxy 
groups that form a link between poorly defined "field" envi- 
ronments and X-ray luminous/optically-rich clusters. If, as 
anticipated, X-ray clusters occupying this luminosity range 
display X-ray temperatures T < 4 keV, they are more likely 
to display the effects of non-gravitational energy input into 
the Intra-Cluster Medium (ICM) than hotter, more mas- 
sive clusters (Ponman, Cannon & Navarro 1999) . Deviations 
of X-ray scaling relations from simple, self-similar expec- 
tations have been studied for structures displaying a rela- 
tively wide range of mass/temperature scales at z < 0.2 
(e.g. Sanderson et al. 2003). However, the study of X-ray 
emitting structures - selected over an extended tempera- 
ture range - at z > 0.2 will provide an important insight 
into the evolution of their X-ray emitting gas. Though de- 
tailed X-ray studies of galaxy clusters at 0.2 < z < 0.6 are 
in progress (e.g. Lumb et al. 2004), few examples of X-ray 
emitting structures displaying temperatures T < 3 keV arc 
currently known at such redshifts. Clearly, compilation of 
a sample of X-ray emitting galaxy groups and clusters to 
z < 0.6 will greatly increase the range of X-ray gas temper- 
atures over which evolutionary effects can be studied. 

Low mass clusters and groups are predicted to be sites of 
continuing galaxy evolution at z < 1 (e.g. Kaufmann 1996; 
Baugh et al. 1996). When identifying such systems, it is im- 
portant to note that extended X-ray emission arises from 
gravitationally bound structures. This is an important dif- 
ference when X-ray selected cluster samples are compared 
to optical/NIR selected cluster samples - whose dynamical 
state can only be assessed with additional velocity data. In 
addition, the X-ray properties of galaxy structures (lumi- 
nosity and temperature) constrain the gravitational mass of 
the emitting structure. Extending the currently known sam- 



ple of galaxy groups and low-mass clusters at z > 0.2 via 
X-ray observations will provide an important group/cluster 
sample with consistent mass ordering. A mass ordered clus- 
ter sample will permit several detailed studies of galaxy evo- 
lution at look back times > 3 — 4 Gyr to be undertaken; 
e.g. morphological segregation and merger-related effects 
(Heldson and Ponman 2003), Butcher- Oemler effects (An- 
dreon and Ettori 1999) and the evolution of colour and lu- 
minosity functions (Andreon et al. 2004a). To date, detailed 
galaxy evolution studies of moderately distant z > 0.2, 
X-ray selected clusters have been performed typically for 
only the most X-ray bright (i.e. massive) systems, e.g. 
L x ([0.3 - 3.5]keV) > 4 x 10 44 ergs s -1 (Yee et al. 1996; 
Dressier 1999). Clearly an improved sample of systems cov- 
ering an extended mass interval will permit a detailed inves- 
tigation of galaxy evolution effects as a function of changing 
environment. 



The X-ray Multi-Mirror (XMM) Large Scale Struc- 
ture (LSS) survey (Pierre et al. 2004) is a wide area X- 
ray survey with the XMM facility with the primary aim 
to extend detailed studies of the X-ray cluster correlation 
function, currently determined at z < 0.2 as part of the 
REFLEX survey (Schuecker et al. 2001), to a redshift of 
unity. However, the XMM-LSS survey features a number of 
secondary aims including determination of the cosmologi- 
cal mass function to faint X-ray luminosities, the evolution 
of cluster galaxy populations and the evolution of the X- 
ray emitting gas in clusters selected over a range of mass 
scales. The nominal point source flux limit of XMM-LSS 
is 5 x 10 -15 ergs -1 cm -2 in the [0.5-2] keV energy band. 
Refregier, Valtchanov and Pierre (2002) demonstrate that 
(assuming a reasonable distribution of cluster surface bright- 
ness profiles) the approximate flux limit for typical extended 
sources is 8 x 10~ 15 ergs -1 cm -2 which corresponds to a X— 
ray luminosity Lx = 1-2 x 10 43 ergs s -1 for a cluster lo- 
cated at a redshift z = 0.6 and Lx = 4.2 x 10 43 ergs s -1 
for a cluster located at a redshift 2 = 1 within the adopted 
cosmological model (see below). 

The above aims are predicated upon the compilation 
of a large, well-defined cluster catalogue. This paper de- 
scribes the first results of the XMM-LSS survey at iden- 
tifying X-ray emitting clusters at z < 0.6. The first clus- 
ters identified at z > 0.6 are presented in Valtchanov et al. 
(2004). The current paper is organised as follows; Section 
2 summarises the X-ray and optical imaging data and the 
methods employed to select candidate clusters. Section 3 de- 
scribes spectroscopic observations and reductions performed 
for a subset of candidate clusters. Section 4 presents the de- 
termination of cluster spectroscopic properties (redshift and 
line-of-sight velocity dispersion). Section 5 presents the de- 
termination of confirmed cluster X-ray properties (surface 
brightness and temperature fitting). Section 6 presents the 
current conclusions for the properties of the initial z < 0.6 
sample. Throughout this paper a Friedmann-Robertson- 
Walker cosmological model, characterised by the present- 
day parameters !)m = 0.3, Qa = 0.7 and Ho = 70 kms -1 
Mpc -1 , is assumed. Where used, h is defined ash = -ffo/(100 
kms -1 Mpc -1 ). 
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2 IDENTIFYING X-RAY CLUSTER 
CANDIDATES 

2.1 X-ray data reduction and source detection 

Galaxy cluster targets presented in this paper were selected 
from a mosaic of overlapping XMM pointings covering a 
total area of approximately 3.5 deg 2 . This data set repre- 
sents all XMM-LSS pointings received by August 2002 and 
includes 15 A0-1 10 ks exposures and 15 Guaranteed Time 
(GT) 20 ks exposures obtained as part of the XMM Medium 
Deep Survey (MDS). 

XMM observations were reduced employing the XMM 
Science Analysis System (SAS) tasks emchain and epchain 
for the MOS and pn detectors respectively. High background 
periods induced by soft-proton flaring were excluded from 
the event lists and raw photon images as a function of energy 
band were created. The raw images for each detector were 
processed employing an iterative wavelet technique and a 
Poissonian noise model with a threshold of 10~ 3 (equivalent 
to 3cr for the Gaussian case) applied to select the significant 
wavelet coefficients (Starck & Pierre 1998). Each wavelet 
filtered image was exposure corrected and an image mask 
(including deviant pixels, detector gaps and non-exposed 
detector regions) was created. 

Source detection was performed on the wavelet fil- 
tered X-ray images em ploying the SExtractor package 
l|Bertin fc Arnouts 19 96'). The discrimination between ex- 
tended (cluster) and point-like sources (mostly AGN) was 
achieved employing a two-constraint test based on the half- 
energy radius and the SExtractor stellarity index of the 
sources. The applied procedure is the optimum method 
give n the XMM PSF and the Poisson nature of the sig- 
nal l|Valtchanov et al. 200lF) . The measurement of extended 
source properties was performed on the EPIC/pn images as 
they provide the greatest sensitivity. The EPIC/MOS im- 
ages were used to discard possible artefacts resulting from 
edge effects associated with the pn CCDs. 

2.2 Optical imaging 

The selection of potential galaxy members within each can- 
didate cluster was performed employing moderately deep 
BVRI images from the CFH12k camera on the Canada 
France Hawaii Telescope (CFHT) obtained as part of the 
VIRMOS deep imaging survey (McCracken et al. 2003) . Ob- 
servations were processed employing the Terapix data re- 
duction pipeline to produce an astrometric and photomet- 
ric image data set. Object catalogues were produced using 
SExtractor. Catalogue detection thresholds as a function of 
photometric band are displayed in Table Q 

2.3 Selecting candidate clusters and member 
galaxies 

The identification of galaxy clusters over an extended red- 
shift interval in X-ray images is limited by the ability of the 
XMM facility to identify extended cluster X-ray emission 
in a lOks exposure. The Half-Energy Width (HEW) of a 

1 htp: //terapix. iap.fr 



Table 1. Detection thresholds as a function of photometric pass- 
band for objects detected in the optical images. 



Filter 


Detection threshold 




(50% AB magnitude completeness 




limit for stellar sources) 


B 


26.5 


V 


26.0 


R 


26.0 


1 


25.4 



on-axis point source is approximately 15" at 1.5 keV. How- 
ever, the HEW displays marked local variations resulting 
from off-axis angle, vignetting and detector gaps. The on- 
axis HEW corresponds to a projected transverse distance of 
120 kpc at a redshift z = 1 within the assumed cosmological 
model. Although this HEW is sufficient to resolve the ex- 
tended emission from massive galaxy clusters to z>l, the ef- 
fect of low central surface brightness, leading to a truncated 
detectable cluster extension, can lead to a cluster being er- 
roneously identified as an unresolved object. Therefore both 
distant clusters and intrinsically compact clusters at all red- 
shifts may potentially appear as only marginally resolved or 
unresolved sources in XMM-LSS X-ray mosaics. A quanti- 
tative assessment of the cluster X-ray selection function will 
form the subject of a future paper (Pacaud et al. 2005). 

The XMM-LSS incorporates a number of different solu- 
tions to the problem of cluster identification, e.g. correlation 
of extended X-ray sources with optical galaxy structures 
(this paper and Valtchanov et al. 2004), investigation of the 
X-ray properties of optically identified structures (and vice 
versa; Andreon et al. 2004b) and the investigation of ex- 
tended X-ray sources lacking optical counterparts together 
with unresolved X-ray sources associated with faint optical 
structures (Andreon et al. 2004b). 

Analysis of the first 3.5 deg 2 of the XMM-LSS survey 
led to the identification of 55 extended sources with fluxes 
greater than 8 x 10 -15 ergs -1 cm -2 which are extended ac- 
cording to the criteria detailed in Section 12.11 The optical 
imaging data corresponding to a 7' x 7' field 2 centred upon 
the location of each extended X-ray source was analysed 
for the presence of a significant galaxy structure showing 
a well defined, red colour sequence (Stanford et al. 1998; 
Kodama et al. 1998). Candidate cluster members were se- 
lected by inspection of the available BVRI photometry of 
objects identified within the field of each extended X-ray 
source. Colour magnitude thresholds were applied interac- 
tively in order to enhance galaxy structures when viewed 
in a VRI pseudo-colour image of each field with X-ray 
contours superposed from the wavelet-filtered X-ray im- 
age. One exception to this procedure was cluster candidate 
XLSS J022722. 3-032141 (see Tabled; the optical data for 
this candidate cluster consisted of the /-band VLT/FORS2 

2 A practical limit determined by the field of view of the multi- 
object spectroscopic facilities employed to observe cluster candi- 
dates (see Section|2J. The field size is large compared to the extent 
of the X-ray emission and obviates any requirement to adjust the 
field centre to maximise the number of candidate cluster mem- 
bers. 
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pre-image obtained to define slit locations. Candidate clus- 
ter members were selected to include all galaxies up to 1.5 
magnitudes fainter than the bright (/ = 17.1), central galaxy 
associated with the extended X-ray source. 

The sample of candidate cluster members generated by 
the above procedures was used to design spectroscopic masks 
for each candidate cluster field. Two multi-slit masks were 
created for each candidate cluster with brighter galaxies 
given higher priority in the slit assignment procedure. Un- 
used regions of each multi-slit mask were employed to sam- 
ple the population of unresolved X-ray sources with bright 
(R < 23) optical counterparts. Further discussion of this 
additional sample will appear elsewhere. 



3 SPECTROSCOPIC OBSERVATIONS 

The candidate cluster sample was observed by the Las Cam- 
panas Observatory Baade telescope with the Low Dispersion 
Survey Spectrograph (LDSS2) during 4-5 October 2002 and 
the European Southern Observatory Very Large Telescope 
(VLT) with Focal Reducing Spectrograph (FORS2) during 
9-12 October 2002. Each instrument (LDSS2 and FORS2) 
is a focal reducing spectrograph with both an imaging and 
a multi-object spectroscopy (MOS) capability. In each case, 
MOS observations are performed using slit masks mounted 
in the instrument focal plane. Details of which cluster candi- 
date was observed with which telescope plus instrument con- 
figuration are provided in Table H The effect ive wavelength 
interval, pixel sampling and spectral resolution generated by 
each instrument combination are indicated in Table [S] 

Spectroscopic observations were reduced employing 
standard data reduction procedures within IRAF 3 : a zero 
level, flat-field and cosmic ray correction was applied to all 
MOS observations prior to the identification, sky subtrac- 
tion and extraction of individual spectral traces employing 
the apextract package. The dispersion solution for each ex- 
tracted spectrum was determined employing HeNeAr lamp 
exposures and all data spectra were resampled to a linear 
wavelength scale. A single spectrophotometric standard star 
from the atlas of Hamuy et al. (1992, 1994) was observed 
during each night and was employed to correct for the rel- 
ative instrumental efficiency as a function of wavelength. 
Removal of the relative instrumental efficiency as a func- 
tion of wavelength does not affect the later determination 
of galaxy redshifts via cross-correlation analysis. However, 
it does permit spectra to be displayed on a relative spec- 
tral flux scale that aids the visual assessment of low quality 
spectra. 

3.1 Spectral classification and redshift 
determination 

In order to confirm the redshift of each candidate cluster, 
the spectroscopic sample generated for each field was con- 
structed to maximise the number of potential cluster mem- 

3 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



bers according to the photometric criteria described in Sec- 
tion [53] Though these criteria were constructed in order to 
identify the characteristic colour signature generated by an 
overdensity of early-type galaxies at a particular redshift, it 
is probable that the spectral sample generated for each can- 
didate cluster is contaminated by the presence of galaxies 
within the target field that are (gravitationally) unassoci- 
ated with the cluster and with stars misidentified as galaxies. 
In order to address this issue and to classify each candidate 
cluster member, all reduced spectra were inspected visually 
to identify contaminating stars and to provide an initial 
estimate of galaxy redshifts based upon the identification 
of prominent features. Individual spectra were then cross- 
correlated with a representative early-type galaxy template 
(Kinney et al. 1996) employing the IRAF routine xcsao 
(Tonry & Davis 1979) . The cross-correlation procedure was 
performed interactively in order to improve the identifica- 
tion of a reliable cross-correlation peak. Spectral regions 
corresponding to the observed locations of prominent night 
sky emission features and regions of strong atmospheric ab- 
sorption were masked within the cross-correlation analysis. 
Computed redshift values have not been corrected to a he- 
liocentric velocity scale. 

Errors in the cross-correlation velocity returned by 
xcsao are computed based upon the fitted peak height 
and the antisymmetric noise component associated with 
the identified cross-correlation peak (Tonry & Davis 1979; 
Heavens 1993). The typical median velocity error computed 
for spectra observed with each spectrograph combination 
described in Table is 75 kms -1 for FORS2+600RI and 
150 kms" 1 for FORS2+300V and LDSS2/medium-red. The 
random error associated with uncertainties in the disper- 
sion solution applied to each spectrum was characterised 
via determination of the error in the wavelength location of 
prominent emission features in the night sky spectrum as- 
sociated with each data spectrum when compared to their 
reference values. The distribution of wavelength residuals 
are considered in velocity space and the rostat statistics 
package (Beers, Flynn & Gebhardt 1990; see Section is 
employed to calculate the bi-weight mean and scale for each 
cluster field. Typical values (for all spectrograph settings) 
of the mean wavelength shift and dispersion computed via 
this method are < ±100 kms -1 and < 50 kms -1 respec- 
tively. Radial velocities (from which redshift of each galaxy 
is determined) are corrected for the mean velocity residual 
in each field and we further assume that the distribution of 
errors in the cross-correlation velocities and in the disper- 
sion solutions to be Gaussian and therefore calculate a total 
uncertainty in the corrected radial velocity by combining 
these two sources of error in quadrature. 



4 DETERMINATION OF CLUSTER 
SPECTROSCOPIC PROPERTIES 

The nature of each candidate cluster contained within the 
spectroscopic sample is assessed employing the available X- 
ray and optical images and the spectroscopic information 
accumulated for each field. The field of each candidate clus- 
ter is inspected visually employing a composite image con- 
taining the CFH12k J?-band greyscale image, X-ray con- 
tours derived from the wavelet-filtered XMM mosaic and 
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Table 2. Observing log of candidate groups and clusters. 



ID a Cluster Right Ascension'' Declination Instrument Grism + Filter # of masks Total exposure time 

per mask (seconds) 



006 


XLSSUJ022145. 


.2- 


-034614 


02: 


21 


:45. 


22 


-03 


:46 


:14. 


1 


FORS2 


300V 4 


- GG435 


2 


2700- 


-2400 


007 


XLSSUJ022406. 


.0- 


-035511 


02: 


24 


:05. 


95 


-03 


:55 


:11. 


.4 


FORS2 


600RI - 


h GG435 


2 


2400- 


-2400 


008 


XLSSUJ022520. 


.7- 


-034800 


02: 


25 


:20. 


71 


-03 


:48 


:00. 





FORS2 


600RI - 


h GG435 


2 


1200- 


- 1200 


009 


XLSSUJ022644. 


.2- 


-034042 


02: 


26 


:14. 


21 


-03 


:40 


:41. 


.8 


FORS2 


300V -| 


- GG435 


2 


1800- 


- 1800 


010 


XLSSUJ022722. 


.2- 


-032137 


02: 


27 


:22. 


16 


-03 


:21 


:37. 





FORS2 


600RI - 


f GG435 


1 


600 


012 


XLSSUJ022827. 


.5- 


-042554 


02: 


28 


:27. 


17 


-04 


:25 


:54. 


3 


LDSS2 


medium-red 


2 


1800- 


- 1800 


013 


XLSSUJ022726. 


.0- 


-043213 


02: 


27 


:25. 


98 


-04 


:32 


:13. 


1 


LDSS2 


medium-red 


2 


900- 


-900 


014 


XLSSUJ022633. 


.9- 


-040348 


02: 


26 


:33 


.87 


-04 


:03 


:48. 





LDSS2 


medium-red 


2 


1800- 


- 1800 


016 


XLSSUJ022829. 


.0- 


-045932 


02: 


28 


:29. 


03 


-04 


:59 


:32. 


2 


LDSS2 


medium-red 


2 


600- 


-900 


017 


XLSSUJ022628. 


.2- 


-045948 


02: 


26 


:28. 


19 


-04 


:59 


:48. 


1 


LDSS2 


medium-red 


2 


1800- 


- 1800 


018 


XLSSUJ022401. 


.5- 


-050525 


02: 


24 


:01. 


16 


-05 


:05 


:24. 


8 


LDSS2 


medium-red 


2 


2400- 


-2400 


020 


XLSSUJ022627. 


.0- 


-050008 


02: 


26 


:27. 


08 


-05 


:00 


:08. 


4 


LDSS2 


medium-red 


2 


1800- 


- 1800 



a In the following text, all clusters are referred to via the reference XLSSC plus the identification number, e.g. XLSSC 006, etc. Clusters 
001-005 correspond to redshift z > 0.6 clusters presented in Valtchanov et al. (2004). 
6 Positions are J2000.0. 



Table 3. Instrumental characteristics for each spectrograph configuration employed during the observations 



Instrument 


Grism 


-I- Filter 


Wavelength 


Pixel 


sampling 


Spectral 








interval (A) 


(A 


pix" 1 ) 


resolution 0, (A) 


FORS2 


300V + 


GG435 


4000- 


-9000 




3.2 


14 


FORS2 


600RI 4 


- GG435 


5000- 


-8500 




1.6 


7 


LDSS2 


medium-red 


4000- 


-9000 




5.1 


14 



a Estimated for each spectrograph via the mean full— width at half-maximum of the HeI5876 arc emission line. All spectral observations 
were performed with a slit width of approximately 1'.'4. 



the available redshifts of all galaxies contained within the 
field (see Figure 0. The redshift distribution generated by 
all spectroscopic redshifts obtained in each candidate cluster 
field is also displayed. To illustrate redshift space overdensi- 
ties, the redshift density computed by applying an adaptive 
kernel (Silverman 1986) to the redshift data is also shown. 
This process provides an initial estimate of the cluster red- 
shift via the identification of three dimensional (position 
and redshift) structures associated with the extended X-ray 
emission. This redshift estimate is then employed to select 
the corresponding peak in the redshift histogram of the field. 
In the case of clusters XLSSC 017 and 020, the redshift of 
each extended X-ray source was determined by determin- 
ing the spatial barycentre of each redshift peak displayed in 
Figure 1 and assigning the redshift grouping closest to each 
X-ray source as the cluster redshift. 

The typically small number (< 20) of objects observed 
spectroscopically in each cluster field limits the usefulness of 
any assessment of spectroscopic completeness and redshift 
confirmation frequency. However, the spectroscopic redshift 
reported for each cluster is in excellent agreement (Az < 
0.02) with the redshift determined independently from the 
location of the red envelope of the corresponding cluster 
colour magnitude relation in R — z colour space (Andreon 
et al. 2004a). 

The sample of cluster members is selected in radial ve- 



locity space employing an iterative method similar to that 
of Lubin, Oke & Postman (2002): the initial cluster sample 
is selected to lie within the redshift interval Az = ±0.06 
of the estimated cluster redshift. Radial velocities relative 
to the cluster centre are calculated within the cluster rest 
frame, i.e. Av = c(z — z)/(l + z) where z is the median 
redshift within the specified interval. The bi-weight mean 
and scale of the radial velocity distribution within this in- 
terval is computed using rostat and galaxies that display a 
velocity difference relative to the central location of greater 
than 3500kms -1 or three scale measures, are rejected and 
the statistical measures recalculated. This procedure is re- 
peated until no further galaxies are rejected. Errors in the 
bi-weight mean and scale are estimated employing a boot- 
strap or jacknife calculation with 10,000 resamplings for 
clusters with greater than or less than 10 confirmed mem- 
bers respectively. Estimates of the bi-weighted mean radial 
velocity and line-of-sight velocity dispersion of each clus- 
ter, and the associated uncertainties, are corrected for biases 
arising from measurement errors employing the prescription 
of Danese et al. (1980). 

Although galaxy redshift and positional information is 
compared to the location of the X-ray source to define the 
initial cluster redshift input to the velocity search algorithm, 
no additional spatial filtering of potential galaxy cluster 
members is performed. However, the projected transverse 
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Figure 1. Optical images, X-ray contours and spectroscopic information for each cluster presented in Table |3] Left panel: adaptively 
smoothed density profile versus rcdshift. Individual rcdshifts are marked by short vertical lines below the density curve. A short vertical 
line above the density curve marks the redshift location of the cluster. The inset histogram displays the rest-frame velocity distribution 
of confirmed cluster members. Right panel: a 7' X 7' ij-band image centred on the extended X-ray source. Wavelet filtered X-ray contours 
are overplotted and squares indicate the position of confirmed cluster members. The X-ray contours typically run from two times the 
background level in each frame to 5 photons/pixel with 10 logarithmic levels. Note that the optical images form a heterogeneous data 
set and are presented to indicate the visual appearance of each cluster. All optical images are orientated with North up and East left. 
Where more than one cluster is detected in the same field a black arrow indicates the location of cluster being considered. 
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Figure 1. Continued. 
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Figure 1. Continued. 
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Table 4. Spectroscopic properties of all redshift z < 
0.6 groups and clusters 



Cluster Redshift a # of members cr v 



(kms 1 ) 


XLSSC 006 


0.429 


39 


821+Jl 


XLSSC 007 


0.558 


10 


oqo + 178 
ozo -191 


XLSSC 008 


0.298 


11 


351±|g 


XLSSC 009 


0.327 


13 


232+|° 


XLSSC 010 


0.329 


8 


420 ± 72 


XLSSC 012 


0.433 


13 


694^04 


XLSSC 013 c 


0.307 


5 


N/A 


XLSSC 014 


0.344 


8 


416 ± 246 


XLSSC 016 


0.332 


5 


703 ± 266 


XLSSC 017 


0.382 


7 


571 ± 282 


XLSSC 018 


0.322 


12 


342+^ 4 


XLSSC 020 


0.494 


6 


265+iM 



a The uncertainty associated with cluster redshifts is 
less than Az = 0.001 in all cases. 

" Velocity dispersion uncertainties are quoted at the 
68% confidence level. 

c The available data for cluster 013 does not generate 
a well-defined velocity dispersion. 

distances sampled by the detector fields of view employed 
to perform the spectroscopy vary between 1.5 and 2 h -1 Mpc 
for clusters located over the redshift interval 0.3 < z < 0.6. 
Studies of velocity dispersion gradients in both local (Girardi 
et al. 1996) and distant (Borgani et al. 1999) X-ray clusters 
- albeit hotter/more luminous systems than those presented 
in the current paper - indicate that integrated velocity dis- 
persion profiles typically converge within radii r < 1 — 2h _1 
Mpc of the X-ray cluster centre. Upon initial inspection, 
the cluster regions sampled by the projected field of view 
sampled by each telescope plus spectrograph combination 
would appear to be well matched to the convergent velocity 
profile of typical hot/luminous clusters. However, inspection 
of Figure indicates that cluster galaxies are typically con- 
firmed within the central regions in each field - a strategy 
required by the necessity to confirm the redshift of galaxy 
structures near the X-ray source. Though the XMM-LSS 
clusters presented in this paper are typically cooler (i.e. less 
massive) than those presented by Girardi et al. (1996) and 
Borgani et al. (1999), and may reasonably be assumed to be 
intrinsically less extensive, an unknown and potentially sig- 
nificant uncertainty is associated with the assumption that 
the computed velocity dispersion figures represent the con- 
vergent velocity dispersion for each cluster. The resulting 
spectroscopic properties of all cluster candidates listed in 
Table Hare given in Table H 



5 DETERMINATION OF GROUP AND 
CLUSTER X-RAY PROPERTIES 

In order to determine the nature of the spectroscopically 
confirmed clusters presented in this paper, additional anal- 
yses of the available XMM data were performed to charac- 
terise the spatial and spectral properties of the X-ray emit- 



ting gas. When combined with the optical redshift and line- 
of-sight velocity dispersion (where available) information, 
these X-ray measures permit a comparison of the cluster 
sample with lower redshift samples and in particular permit 
the approximate mass interval occupied by z < 0.6 XMM- 
LSS clusters to be understood. 



5.1 Morphological properties 

The X-ray surface brightness distribution of each cluster 
was modelled employing a circular /3-model, of the form 



where the coordinate r is measured in arcseconds with re- 
spect to the centre of the X-ray photon distribution, ro is 
the core radius, A is the amplitude at r — 0, and a = 3 13— i. 

Images and exposure maps for each cluster field were 
created for the three EPIC instruments (MOS1, MOS2 and 
pn) separately in the [0.5-2] keV energy band. Images of the 
appropriate PSF were created (using S AS-calview) , with the 
appropriate energy weighting together with the off-axis and 
azimuthal angles of the source. Square regions (of sizes rang- 
ing between 175" and 500" on a side) were selected around 
each source and around a nearby, source-free, background 
region. Mask images were also created and employed to re- 
move from further analysis regions associated with chip gaps 
and serendipitous point sources lying close to each cluster 
source. 

We used the Sherpa package from the CI AO analysis 
system to fit a model of the form given by Equation 1 to the 
X-ray data for each EPIC instrument. The quality of fit to 
the three instruments was optimised using the Cash statistic, 
providing a maximum likelihood fit, which accounts properly 
for the Poissonian nature of the data. Each model incorpo- 
rates a flat background model (where the background level 
is determined employing the associated background region) 
and a /3-model convolved with the appropriate PSF. For 
each of the three instrument models determined for each 
cluster, the values of the core radius ro, position (xo, yo), 
and slope a were constrained to be identical. Only the nor- 
malisation for each model was permitted to vary. 

Table |K| lists the best-fit structural parameters for each 
of the confirmed cluster sources. Figure shows flux con- 
tours corresponding to the best-fitting surface brightness 
model overplotted on the X-ray emission for each candidate. 
The most interesting outcome from these fits is the low value 
of the fitted /3 parameter for many of these systems, com- 
pared with the typical value of (3 = 0.66 determined for clus- 
ters (Arnaud & Evrard 1999). The median value of (3 = 0.45 
found here agrees with the value of 0.46 derived for a sample 
of low redshift groups by Heldson & Ponman (2000). Tabled 
also lists the value of the normalised \ 2 statistic computed 
from a comparison of the best fitting model to the radi- 
ally averaged surface brightness profile for each source. The 
value of the \ 2 statistic for each cluster generally indicates 
that the computed model provides a statistically acceptable 
fit. The two systems with computed \ 2 values significantly 
less than one (14 and 16) display some of the lowest count 
levels in the sample and, in these cases, the radially aver- 
aged counts may be better described by a Poissonian rather 
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Figure 2. Model surface brightness distributions for the X-ray group and cluster sample. Contours represent the best-fitting model con- 
volved with the detector response and are overplotted on the full (pn+MOSl+MOS2) Gaussian smoothed X— ray emission corresponding 
to the [0.5—2] keV energy band. The Gaussian smoothing scale (i.e. the standard sigma) is 5 pixels or 12'.'5.The linearly spaced contours 
represent the MOSf model (only the normalisation varies between the three instrument models), except for XLSSC 007 and 013, where 
the pn model is displayed. The images are of size 'boxside' (see Table l5l . 



Table 5. Morphological X-ray parameters determined for con- 
firmed groups and clusters. Displayed errors are la and, in the 
case of cluster 19, no reliable error information could be deter- 
mined. Values of x 2 were computed over a sum of radial bins 
for each cluster. These values are provided to indicate the overall 
merit of each fit. 



Cluster 


boxside 


ro 





x 2 




arcsec 


arcsec 




(per d.o.f.) 



XLSSC 006 
XLSSC 007 
XLSSC 008 
XLSSC 009 
XLSSC 010 
XLSSC 012 
XLSSC 013 
XLSSC 014 
XLSSC 016 
XLSSC 017 
XLSSC 018 
XLSSC 020 



250 
400 
250 
250 
400 
500 
175 
175 
200 
200 
250 
200 



24 0+ 14 - 
^+8.7 
-6.0 
1+52.0 
-19.0 
;+3.3 
-2.6 
-,+10.0 

-7.2 
S +14.9 

-6.2 
7+39 
-4.7 
^+3.5 
-2.4 

i2.i±L 6 

3+8.1 



9.6: 
44.0;! 

8.5j 
29.2^ 
11.1 

4.7! 

2^ 



h.04 
-.02 
1+.07 
-.06 
K09 
.06 

0.67+' 73 



0.58Z 
0.40^ 
0.44 4 



0.44" 



.17 
.03 
.()_> 
.08 



24. 8^ 



0-54_. 48 
0.67±% 
0A0±f 6 

- + .10 
.07 

0.55±;°| 

42+ 03 
u -^ z -.03 

h.04 



0.45^ 



22.2 



.16 



1.34 
1.15 
1.26 
0.96 
1.68 
1.39 
1.34 
0.57 
0.79 
1.11 
1.14 
0.99 



than Gaussian noise distribution - partially invalidating the 
application of a \ 2 merit function. 



5.2 Spectral properties 

X-ray spectra for each cluster were extracted within a 
source circle of radius ranging between 30" and 90". The 
corresponding background spectrum was extracted from a 
surrounding annulus. Sources adjacent to the cluster were 
nagged and removed from the spectral analysis employing 
the source region file generated by the original source ex- 
traction software. Source mapping used the stacked pn + 
MOS1 + MOS2 image of each cluster field in the [0.5-2] keV 
band only. The source extraction and background regions 
applied to XLSSC 013 is shown in Figure [2] as an example 
of the procedure. 

Extracted spectral data corresponding to pn + MOS1 
+ MOS2 detectors were fitted simultaneously. The fitting 
model consists of an absorbed APEC hot plasma model 
(Smith et al. 2001) with a metal abundance ratio set to 
Grevesse and Sauval (1999) values. The hydrogen absorp- 
tion is modeled using a wabs model with Nh fixed at the 
Galactic value, i.e. Nh ~ 2.6 x 10 20 cm" 2 . Spectral data 
are resampled such that the associated background spec- 
trum displays 5 counts per bin and model values are com- 
pared to the data by computing the corresponding value of 
the C-statistic (see Appendix lAl for a justification for this 
approach). Model fitting is performed in two stages; first 
the temperature and abundance are fixed (T = 0.5 keV, 
Z/Zq = 0.3) and only the count normalization is fitted. 
Once a best-fitting spectrum normalisation has been com- 
puted, the best-fitting temperature is computed assuming a 
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Figure 4. Spectral data for groups and clusters described in Table [B] for which a temperature was fitted. For each system the upper 
panel depicts the spectral photon flux for the MOS1 (thick black crosses), MOS2 (thin black crosses) and pn (grey crosses). The data 
have been resampled to 20 photons per spectral bin for display purposes only. The spectral model applicable to each detector is also 
shown by the appropriately coloured solid line, i.e. MOS1 (solid black line), MOS2 (dotted black line) and pn (grey dashed line). The 
lower panel shows the deviation of the data from the model for each detector expressed in units of normalised x 2 per spectral bin. 
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Figure 3. Source and background extraction regions applied to 
XLSSC 013. The greyscale shows the photon map generated by 
the stacked pn + MOS1 + MOS2 image. The central solid circle 
represents the source extraction region. Solid circles crossed with 
a diagonal line indicate regions excluded from the source aperture. 
The two large concentric dashed circles define the background 
annulus and dashed circles crossed with diagonal lines indicate 
regions excluded from the background aperture. 

fixed metal abundance. The fitting results are displayed in 
Table |S| Table |S| also shows the value of the normalised C- 
statistic computed from a comparison of the spectral data 
for each cluster to the quoted model. In each case a statis- 
tically reasonable agreement is obtained. Spectral data for 
all clusters possessing a fitted temperature in Table |f>] arc 
displayed in Figure 2] In the cases of clusters 7, 14, 16, 17 
and 20, no spectral fit was possible and a temperature of 1.5 
keV (typical of the sample as a whole) was assumed for the 
purpose of calculating the flux from the system. 

X-ray luminosities for each cluster source have been 
determined within a uniform physical scale derived from the 
cluster overdensity radius. In the present study we employ 
the radius rsoo, within which the total mean density of the 
system is 500 times the critical density of the universe at 
the redshift of the system. The value of rsoo is computed 
using an isothermal /3-model (Ettori 2000) and employing 
the fitted gas temperature and /3 value for each system 4 . 

We have compared the value of rsoo derived from the 
above isothermal model for each cluster to the correspond- 
ing values obtained using the observed relationship between 
cluster mass and temperature. Employing the data pre- 
sented by Finguenov, Reiprich and Bohringer (2001), based 
on systems ranging from Tx = 0.75 — 14 keV, we obtain the 
following relation 

r 500 = 0.391 t£ 63 /i 7 o(^) _1 Mpc, (2) 

4 We modify Equation A2 of Ettori (2000) to reflect the varying 
redshift dependence of Hubble parameter in a matter plus lambda 
cosmological model. 



where hjo(z) describes the redshift evolution of the Hubble 
parameter (scaled to 70 kms -1 Mpc -1 ) in the assumed cos- 
mological model. The value for rsoo for each cluster derived 
using the above methods typically agree to within ±10% 
with the maximal difference being ±20%. Based upon this 
comparison we employ rsoo values based upon the isother- 
mal model in the rest of this paper. 

In those cases where a successful spectral fit was ob- 
tained, the derived bolometric source flux was extrapolated 
to rsoo using the /3-model determined in Section 5.1. For 
the four systems with detected flux but with no success- 
ful spectral fit, the corresponding best-fitting spatial model 
was employed to compute the bolometric flux within rsoo, 
assuming a T = 1.5 keV emission spectrum. Uncertainties 
on the resulting value of bolometric luminosity are available 
only for those systems with spectral fits, in which cases the 
uncertainty is simply derived from the error on the fitted 
spectral model normalisation. Computation of luminosity 
uncertainties employing this method does not include the 
effects of the uncertain extrapolation of the surface bright- 
ness model to rsoo- The aperture correction factor for each 
cluster, A, defined as the relative change in the integrated 
/3-model profile obtained by varying the integration limit 
from r spec to rsoo, is shown in Tabic H2 



6 THE NATURE OF XMM-LSS SURVEY 
SELECTED GROUPS AND CLUSTERS AT 

Z < 0.6 

In this section we compare the properties of the XMM- 
LSS groups and clusters at z < 0.6 with X-ray group and 
cluster samples in the literature. Figure compares the 
Lx versus Tx distribution of XMM LSS clusters confirmed 
at z < 0.6 to both the distribution formed by the Group 
Evolution Multi-wavelength Study (GEMS) sample of lo- 
cal (z < 0.03) X-ray emitting galaxy groups (Osmond and 
Ponman 2004; hereafter OP04) and the sample of Marke- 
vitch (1998; hereafter M98) containing clusters at z < 0.09 
with ASCA temperatures and ROS AT luminosities 5 . Figure 
[S] indicates that XMM-LSS clusters occupy a region of the 
Lx versus Tx plane ranging from cool (Tx > 0.9 keV), low 
luminosity (Lx(rsoo) > 4 x 10 42 ergs s _1 ) X-ray groups to 
moderate temperature (Tx = 5 keV), moderate luminosity 
(Lx (rsoo) = 4 x 10 44 ergs s _1 ) clusters. Though the sample 
of X-ray systems presented in this paper is not statistically 
complete, it is representative of the properties of the com- 
plete flux-limited sample currently under construction. Due 
to the steeply rising nature of the XLF toward faint X-ray 
systems, it is anticipated that the larger, statistically com- 
plete sample of X-ray structures identified by XMM-LSS at 
z < 0.6 will be dominated by such galaxy group and low 
mass cluster systems. 

It can be seen from Figure^] that our XMM-LSS groups 
and clusters appear to be in good agreement with a linear 
fit to the Lx versus Tx distribution of lower redshift group 

5 Note that, although luminosities in M98 are quoted within 1 
h"loa Mpc apertures and not rsoo as used in this paper, the appro- 
priate correction factors are typically considerably smaller than 
the 5% calibration uncertainty associated with the luminosities 
themselves. 
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Table 6. Spectral X-ray parameters determined for confirmed groups and clusters. Values for exposure time, t exp , and total counts are 
summed over all three detectors. Where the letter "F" follows a tabulated temperature (T) value, this indicates that the value was fixed 
in the fitting procedure. The definition of the aperture correction factor, A, is provided in the text. Displayed errors are lcr. 
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Figure 5. Distribution of X— ray luminosity computed within a 
scale radius rsoo and temperature for all XMM— LSS groups and 
clusters currently identified at 2 ^ 0.6 (solid squares). Also indi- 
cated are values of X-ray luminosity and temperature determined 
for the low redshift group sample of OP04 (open squares) and for 
the cluster sample of Markevitch (1998) (open triangles). The 
solid line indicates an orthogonal regression fit to the Lx versus 
Tx relation for both the group and cluster sample incorporating a 
treatment of the selection effects present in each sample (Heldson 
and Ponman 2005) - see text for details. 



and cluster samples. The local fit to the Lx versus Tx dis- 
tribution takes the form \ogLx = 2.91 logTx + 42.54 and 
is derived from an orthogonal regression fit to the combined 
OP04 plus M98 samples incorporating a treatment of the se- 
lection effects present in each sample (Heldson and Ponman 
2005). The comparison of XMM-LSS groups and clusters to 
this local relationship may be quantified (Figure by cal- 
culating a luminosity enhancement factor, F = L bs/L pre d, 
where L a t, s is the observed cluster X-ray luminosity within 



a radius, rsoo, and L pre d is the luminosity expected applying 
the fitted Lx versus Tx relation computed for the local fit 
and the XMM-LSS measured temperature. Neglecting the 

5 systems assigned a fixed temperature (XLSSC 007, 014, 
016, 017 and 020 - for which the temperature uncertainty is 
unknown) , the median enhancement factor of the remain- 
ing 6 systems is F = 1.09. For comparison, the expected 
enhancement in Lx due to self-similar evolution, scaling to 
rsoo, is a factor of 1.23 at the typical redshift (z = 0.4) of 
our sample. Therefore, given the observed spread in the en- 
hancement values, the observed negative deviation from the 
self-similar expectation is not large. 

Given the modest size and the statistically incomplete 
nature of our current sample, we regard these results as 
in need of confirmation. Ettori et al. (2004) also report 
evolution weaker than the self-similar expectation from a 
sample of 28 clusters at z > 0.4 with gas temperatures 
3 keV < T < 11 keV. The combined effect of self-similar 
scaling with the negative evolution reported by Ettori et al. 
(2004) would result in an enhancement factor F = 0.86 6 at 
a z — 0.4 (see Figure Though the overlap of the Ettori et 
al. (2004) sample and the systems contained in the present 
work is limited, the X-ray luminosities appear to describe 
similar trends. 

The relationship between the specific energy contained 
within the cluster galaxy motions, compared to the X-ray 
emitting gas, is described via the /3 apec parameter 

Papec kTx/[J.m p ^ 

where a v is the line-of-sight galaxy velocity distribution, k 
is the Boltzmann constant, Tx is the X-ray gas temperature 
and fiTrip is the mean particle mass within the gas (Bahcall 

6 Lubin 1994). Figure |7| displays the value of j3 sp ec com- 
puted for the XMM-LSS z < 0.6 cluster sample as a func- 
tion of computed X-ray luminosity extrapolated to a radius 
rsoo (i-e. the consistent measure adopted in this paper). A 

6 This enhancement factor assumes self-similar evolution and an 
additional factor (1 + z) Bz where B z = —1.04, following the 
nomenclature of Ettori et al. (2004). 
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Figure 6. Enhancement factor, F = L a i, a / L pre di computed for 
six XMM-LSS groups and clusters located at 2 ^ 0.6 plotted 
versus the X-ray temperature of each system (see text for ad- 
ditional details). Horizontal lines indicate expected values of F: 
the short dashed line indicates the value F = 1.23 expected from 
self-similar considerations. The dotted line indicates the value of 
F expected at z = 0.4 based upon Ettori et al. (2004). 



number of systems have been excluded from Figure|7] in ad- 
dition to the systems excluded from the computation of the 
luminosity enhancement factor above, XLSSC 013 does not 
possess well defined galaxy velocity dispersion thus prevent- 
ing computation of f3 spe c- Figure |7| indicates typical values 
Pspec < 1 for the XMM-LSS z < 0.6 sample (with the ex- 
clusion of the above mentioned systems), the median value 
of /3 sp ec for this restricted sampled is {(3 ape c = 0.61). This 
value may be compared to values of /3 spec w 1 reported by 
OP04 for luminous X-ray groups (i.e. Lx(rsoo) > 10 42 ergs 

A clear concern when interpreting the trend of low /3 spec 
values is the extent to which the cluster galaxy velocity dis- 
persion estimates may be biased toward lower values. Po- 
tential uncertainties associated with the velocity dispersion 
values presented in this paper are discussed in Section [I] 
However, in the overwhelming majority of clusters observed 
in detail, the integrated velocity dispersion profile of galaxy 
clusters is a decreasing function of projected radius from the 
cluster centre (Girardi et al. 1996; Borgani et al. 1999). If 
the integrated velocity dispersion profiles of the XMM-LSS 
clusters presented in this paper display similar behaviour to 
hotter clusters, then the expectation arising from computa- 
tion of the cluster velocity at some fraction of the convergent 
radius is that the velocity dispersion will be overestimated 
and will result in values of fl spec biased to higher values. The 
extent of any such bias is difficult to quantify in the current 
data set. However, the implication is that the value of f3 spec 
displayed in Figure |7| would not increase with the addition of 
velocity dispersion measurements extending to larger radii. 

The low values of /3 spec apparent in Figurc|7| are similar 
to those seen by OP04 in lower luminosity groups (Lx < 



Figure 7. Values of /3 spec computed for eight XMM-LSS groups 
and clusters at 2 ^ 0.6 (see text for details) plotted versus the 
X— ray luminosity for each system. The horizontal dashed line in- 
dicates the value f3 Bpec = 1. 

10 42 ergs s _1 ) at low redshift. The origin of these low values 
of /3 sp ec is far from clear, but OP04 argue that it appears 
to result primarily from a reduction in a v , rather than an 
enhancement in Tx ■ Whatever the cause, our results provide 
tentative evidence that these effects are operating in hotter 
and more X-ray luminous systems at higher redshift. We 
are currently in the process of conducting magnitude limited 
spectroscopy of a sample of Tx ~ 1 keV systems at z = 0.3 
in order to provide a more robust picture of the dynamics 
of low temperature X-ray systems. 



7 CONCLUSIONS 

We have presented twelve newly identified X-ray selected 
groups and clusters as part of the XMM Large-scale Struc- 
ture Survey. The procedures employed to detect and clas- 
sify sources in X-rays, and to subsequently confirm each 
source via optical imaging and spectroscopic observations 
have been described in detail. 

We have emphasized throughout this paper that the 
current sample of clusters is not complete in any statisti- 
cal sense. The presentation of a larger, complete sample of 
X-ray clusters located at z < 0.6 will form part of a future 
publication. However, the current sample of X-ray clusters 
at z < 0.6 presents a number of interesting features: most 
importantly, the sample is dominated by low X-ray temper- 
ature systems located at redshifts much greater than that 
presented by previous X-ray studies. Such systems are pre- 
dicted to display the effects of pre-heating or additional 
energy input into the ICM to a greater extent than hot- 
ter, more massive systems. The identification of such low- 
temperature systems at look-back times up to 5.7 Gyr pro- 
vides an important baseline over which to study the extent 
to which such systems evolve. 
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We find tentative evidence that these high redshift 
groups are more luminous than local systems, at a given 
temperature, in agreement with recent work on richer clus- 
ters. However, our results suggest that group luminosities 
may be evolving less rapidly than higher temperature clus- 
ters when compared to self-similar models. If this is con- 
firmed to be the case, then the steepening of the Lx-Tx 
relation at low temperatures reported in local samples, may 
continue at higher redshift. We also find preliminary indica- 
tions that the poorly understood tendency for the specific 
energy in the gas to exceed that in the galaxies in poor 
groups, extends to systems with higher values of L x and 
Tx at 2 ~ 0.4. The completion of a larger and statistically 
complete sample of intermediate redshift groups from the 
XMM-LSS survey, should allow these results to be placed 
on a firm statistical footing in the near future. 
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APPENDIX A: FITTING SIMULATED 
CLUSTER SPECTRA 

The current study extends X-ray spectral observations of 
distant galaxy groups and clusters to low integrated signal 
levels (~ 100 photons above the background). It is therefore 
prudent to assess the reliability of temperature measures 
and associated uncertainties computed via model fits to such 
faint spectra by repeating the fitting procedure for a grid of 
simulated spectra created to reproduce the properties of the 
observed sample. 

The source model used to simulate group and clus- 
ter spectra employs an APEC of an optically thin plasma 
(Smith et al. 2001). This model depends upon four param- 
eters: temperature, metal abundance, redshift and a nor- 
malisation representative of the emission integral. To sim- 
ulate galaxy groups and clusters the abundance is set to 
Z/Zq = 0.3 with solar abundance ratios set to Grevesse and 
Sauval (1999) values. The redshift of the simulated source is 
set to z = 0.3, typical of the sources presented in this paper. 
Photoelectric absorption described by a wabs model within 
Xspec using Morrison and McCammon (1983) cross-sections 
was applied to the source model with the neutral hydrogen 
column density fixed to Nh = 2.6 x 10 20 cm -2 - the mean 
value for our sample according to the Hi distribution map 
of Dickey and Lockman (1990). The instrumental response 
was modeled using the redistribution matrix and ancillary 
response files from observations of XLSSC 006. 

A model pn + MOS1 + MOS2 spectra was created 
and a conversion factor applied to generate spectra of the 
required integrated count level over the spectral interval 
[0.3-10] keV in a 10,000 second exposure. Each simulated 
spectrum is generated from this model using Poissonian 
considerations. Each source spectrum is accompanied by 
a background spectrum created from a Poisson realisation 
of a background model of normalisation and shape consis- 
tent with observed cluster backgrounds. Spectra were sim- 
ulated according to this procedure for temperatures equal 
to 1,2,3 and 5 keV. At each temperature, 50 spectra were 
simulated at each point of a grid of integrated count levels 
100, 200, . . . , 1000. 

Spectral fitting follows the same approach as applied 
to observed data, i.e. temperature and spectum normalisa- 
tion are permitted to vary while the abundance is fixed. 



Data from pn + MOS1 + MOS2 are are combined within 
Xspec with the response files from XLSSC 006 and the en- 
ergy range [0.3-10] keV is conserved. The spectral energy 
range corresponding to [7.5-8.5] keV measured by the pn 
detector is ignored as it contains strong instrumental line 
emission that is not well corrected by our data modeling 
process. The best-fitting model is then determind by min- 
imising a modified C-statistic and the 1 a uncertainty about 
the best fit model is computed. While the C-statistic is in- 
tended to work efficiently on unbinned data a comparison 
of the fitted temperature to the input value indicates a ten- 
dency to underestimate the temperature of spectra of in- 
put temperatures < 5 keV displaying count values < 1000 
counts using this procedure. This bias is indicated in Figure 
lAll and appears to arise from the fact that the statistic used 
in Xspec represents a modified C-statistic that accounts for 
statistical fluctuations in the background estimation. Our 
understanding of the problem is that this modified statistic 
fails at estimating model parameters when there is a signif- 
icant number of background bins containing zero photons. 

Resampling the data to prevent the occurence of spec- 
tral bins containing zero counts minimises this negative tem- 
perature bias. The resampling factor is determined by re- 
quiring that the background spectrum associated with each 
source display a specified minimum count level per spectral 
bin. Determining the resampling factor from the background 
spectrum represents a sensible approach as the background 
counts are more numerous and therefore minimise the loss 
of spectral information in the source spectrum. Although 
a small positive bias is introduced to the fitted tempera- 
tures of very low count level spectra (< 300 counts) when 
the data are resampled, the amplitude of this bias is less 
than 10% when the data are resampled to contain 5 counts 
per spectral bin (Figure [All . Applying a larger resampling 
factor increases the positive temperature bias - which can 
be understood in terms of the spectral smoothing that the 
resampling procedure represents. 

We therefore resample the observed data to generate a 
background spectrum containing 5 counts per spectral bin. 
This approach generates fitted temperatures that agree with 
the input temperature to < 10%. In addition, comparison of 
the distribution of fitted temperatures at any given combina- 
tion of temperature and count level to the temperature un- 
certainty returned by Xspec indicates that the Xspec quoted 
errors on fitted temperatures overestimate the distribution 
of fitted temperatures by a factor typically less than 2. Due 
to the various assumptions that enter the simulation proce- 
dure it therefore seems reaonable to provide Xspec quoted 
temperature uncertainties as a conservative error estimate. 

This paper has been typeset from a TjnX/ I^TjhX file prepared 
by the author. 
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Figure Al. A comparison of Xspec computed temperatures for simulated group and cluster spectra employing the C-statistic with 
two different resampling approaches. Left panels indicate the results for unbinned spectra. Right panels indicate the results for binned 
spectra such that the background spectrum displays a minimum of 5 counts per spectral bin. Panels in each row correspond to spectral 
models with the indicated input temperature (also shown by the horizontal dashed line). In each panel, data points represent the mean 
Xspec computed temperature returned from the set of simulated spectra as a function of total input counts. Filled squares plus error 
bars indicate the mean computed temperature and the distribution of temperatures accounting for 68% of the sample. Open squares 
plus error bars indicate the mean computed temperature and the median 1 cr uncertainty returned by Xspec (open squares are shifted 
to the right by 30 counts with respect to the filled squares for clarity). 



